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Abstract 

An asymptotic technique is presented for finding the spectral efficiency of multi-antenna links in wireless 
networks where transmitters have Channel-State-Information (CSI) corresponding to their target receiver. Transmitters 
are assumed to transmit independent data streams on a limited number of channel modes which limits the rank of 
transmit covariance matrices. This technique is applied to spatially distributed networks to derive an approximation 
for the asymptotic spectral efficiency in the interference-limited regime as a function of link-length, interferer density, 
number of antennas per receiver and transmitter, number of transmit streams and path-loss exponent. It is found that 
Qh targeted-receiver CSI, which can be acquired with low overhead in duplex systems with reciprocity, can increase 



(N 



spectral efficiency several fold, particularly when link lengths are large, node density is high or both. Additionally, 
the per-link spectral efficiency is found to be a function of the ratio of node density to the number of receiver 
antennas, and that it can often be improved if nodes transmit using fewer streams. These results are validated for 
finite-sized systems by Monte-Carlo simulation and are asymptotic in the regime where the number of users and 
antennas per receiver approach infinity. 



Index Terms 

O . MIMO, Wireless Networks, Antenna Arrays, Stochastic Geometry, Ad-hoc networks. 

I. Introduction 

Multiple antenna systems are attractive for use in point-to-point and ad-hoc wireless networks since they can 
^ ■ suppress interference, increase data rates by spatial multiplexing or coherent signal combining, and are robust to 
^-j. , channel variations. The increased robustness available through the increased diversity and multiplexing capability 
^ [ of multi-antenna systems in isolated Additive White Gaussian Noise (AWGN) channels has been well studied and 
Q ■ is described in detail in wireless communications texts such as [2]. The major contributions in this area which 
O . addressed the multiplexing capability of multi-antenna systems in AWGN channels include [3], [4], [5]. The trade- 
^7"! ' off between diversity (which increases robustness) and multiplexing (which increases data rates) was studied in 

A central assumption that determines the performance of Multiple-Input-Multiple-Output (MIMO) systems is the 
| availability of Channel-State-Information (CSI) at the transmit(Tx) or receive(Rx) side of links. It is known that 
- - ■ to achieve capacity in Additive -White-Gaussian-Noise (AWGN) channels with Rx CSI and without Tx CSI, the 
transmitter should send equal power, independent data streams on each antenna. With Tx CSI, capacity is achieved 
by parallelizing the channel between the Tx and Rx using a Singular Value Decomposition (SVD) of the channel 
matrix. Most works on multi-antenna systems in the literature assume that receivers have CSI, an assumption which 
is realistic since receivers can estimate channel parameters from signals received from transmitters. 

The performance of multi-antenna systems in the presence of interference has been relatively less well studied 
although the ability of antenna arrays to suppress undesired signals is well known in the signal processing community 
(e.g. see [7]). MIMO communications systems employing interference suppression by transmitters and receivers 
were analyzed under different CSI assumptions by [8] and [9]. 

The main contributions of this article are a technique to find the asymptotic spectral efficiency of multi-antenna 
links in ad-hoc wireless networks with limited CSI at transmitters and arbitrary distribution of interference powers 

S. Govindasamy is with the Franklin W. Olin College of Engineering. D. W. Bliss and D. H. Staelin are with the Research 
Laboratory of Electronics, Massachusetts Institute of Technology (MIT). D. W. Bliss is also with MIT Lincoln Laboratory, (email: 
siddhartan.govindasamy@olin.edu, staelin@mit.edu, bliss@ll.mit.edu). This research was supported in part by the National Science 
Foundation under Grant ANI-0333902. Opinions, interpretations, conclusions, and recommendations are those of the authors and are not 
necessarily endorsed by the United States Government. Portions of this material have been presented at the Asilomar Conference on Signals, 
Systems and Computers 2009, Pacific Grove CA, USA [1]. 



2 



subject to a convergence criteria. In particular, we find an expression for the spectral efficiency in the interference- 
limited regime when interference is due to spatially distributed nodes. Transmitters are assumed to have CSI between 
themselves and their target receivers but to no other receiver in the network. We refer to this form of CSI as Tx-Link 
CSI. We assume that each transmitter is restricted to sending M independent data streams on M channel modes, 
and allow arbitrary correlation between the power allocated to each stream for any given transmitter. The covariance 
matrices of the signals at the antennas of each transmitting node are thus matrices of rank-M. Transmitters are 
assumed to have K > M antennas. Limiting the number of transmit data streams in this fashion is known to 
increase the network spectral efficiency of wireless systems as shown analytically in [10], [11] for systems without 
Tx CSI, and by simulation for systems with Tx CSI in [12]. Note that in some systems, the rank of the transmit 
covariance matrix may be greater than the total number of independent streams, e.g. if spatial repetition coding 
is used. However, in this work, we assume that the number of independent data streams equals the rank of the 
transmit covariance matrices. 

The data rates achievable in wireless networks with interfering MIMO links and no Tx CSI was studied in [10] 
who found that transmitters should transmit a single data stream (i.e. with a rank-1 transmit covariance matrix) from 
one of their antennas in the high interference regime. When interference is low, the links are essentially isolated 
and it is optimal for transmitters to send equal-power streams on each antenna as for AWGN channels. In [13], 
the authors analyzed one-to-one interfering links in the regime where the number of users goes to infinity, both in 
systems without Tx CSI and systems with Tx Link CSI. The authors derived upper bounds to the mean network 
spectral efficiency (b/s/Hz) of such systems, which were found to be constant if the number of receiver antennas per 
node increases linearly with the number of transmitter antennas. Neither [10] or [13] model the spatial distribution 
of nodes however. 

The distribution of nodes in space is key to understanding large wireless networks since signal and interference 
strengths depend on relative node locations. Spatially distributed wireless networks with multiple-antenna links and 
no Tx CSI were studied in [11], which found the asymptotic spectral efficiency as the number of receiver antennas 
./V and interferers in the network n tend to infinity. We found that the mean spectral efficiency was a function of the 
ratio of the number of receiver antennas to the product of node density and the square of the link length, implying 
that constant mean per-link spectral efficiency can be maintained by scaling the number of receiver antennas with 
node density. Recently, [14] and [15] found exact expressions for the Cumulative-Distribution-Function (CDF) 
of the Signal-to-Interference-Ratio (SIR) in spatially distributed networks with MMSE receivers but no CSI at 
the transmitters with [14] considering single-stream transmission and [15] considering multi-stream transmissions. 
Additionally, [16] found that it is possible to scale the network spectral efficiency per unit area linearly with node 
density if the number of receiver antennas scales with node density by using a partial zero-forcing receiver. 

Spatially distributed systems with Tx-CSI have been studied in a number of works such as [17], [18], [19], 
and [20]. In [17] and [18] the receiver did not use its degrees of freedom to mitigate interference whereas [19] 
considered a partial zero-forcing receiver structure. An asymptotic analysis was used in [20] to analyze the spectral 
efficiency of MIMO links in spatially distributed networks of finite area when transmitters use single-streams with 
maximal-ratio-transmission (the transmit-side analogue of MRC), and MMSE receivers. In that work, the authors 
approximate the Signal-to-Interference -plus-Noise-Ratio (SINR) as a gamma distributed random variable and find 
an approximation to the CDF of the SINR as the number of receiver antennas and number of nodes in a finite-area 
network go to infinity. 

In contrast, in addition to Tx-Link CSI and multiple transmit data streams, our work assumes optimal decoding 
at the receiver with interference treated as spatially-colored noise, and provides a general framework which is also 
applicable to systems where interference powers do not depend on node locations, which is relevant to power- 
controlled cellular architectures. Additionally, we assume a constant and finite density of users when applying our 
general technique to spatially distributed networks. 

Our results concerning spatially-distributed networks characterize the spectral efficiency of multi-antenna links 
in ad-hoc wireless networks as a function of tangible parameters such as link length, node density, number of 
antennas, and path-loss exponent, and are useful for system designers to explore the trade offs between increased 
hardware costs of using more antennas or transmit data streams. Additionally, these results enable us to compare 
the spectral efficiency gains that Tx Link CSI provides. 

The asymptotic techniques used here are closely related to several works in the literature such as [21] and 
[22]. In these works, the Signal-to-Interference-plus-Noise-Ratio (SINR) of random Code-Division-Multiple- Access 
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(CDMA) systems, given by terms of the form s'Rs, are shown to converge to asymptotic limits that depend on 
the structure of the limiting distribution of the eigenvalues of the covariance matrix of received interference powers 
R. In both [21] and [22], the vector of signatures s is assumed to have Independent Identically Distributed (IID) 
entries. Here, the SINR associated with a given stream from a transmitter has the form Au^Ru, where u is a 
singular vector associated with a Gaussian random matrix and A is the square of its associated singular value. Since 
u is a unit-norm isotropic vector, its entries are not IID and hence we cannot directly apply the results of [21] 
and [22]. Another related work is [23] which analyzed the limiting SIR of Random CDMA systems with multiple 
antennas. While it allows for some correlation in the entries of the signature vector, its assumptions do not admit 
unit-norm isotropic vectors. Other related works include [24] who analyzed the joint asymptotic SIR distribution of 
multiple transmitters communicating to a single receiver in CDMA systems and [25] who analyze the asymptotic 
capacity in MIMO Multiple-Access and Broadcast Channels with a fixed number of users and number of antennas 
per user going to infinity. 

Furthermore, [21], [22], and [23] analyze the limiting SINR as a function of the limiting distribution of the 
received interference powers, whereas we explicitly consider the distribution of transmit powers and path-losses 
with arbitrary correlation between the power allocated on each stream by the transmitters. While this assumption 
does not exclude the water-filling power allocation, we will explore practical constant-power approaches in the 
applications of the techniques we develop. 

The remainder of the paper is organized as follows. Section [II] presents the network model starting with a general 
model of a network with one-to-one links, and followed by the spatially distributed network model. Section [III] 
contains the main results for the general network model. Section[lV] applies the main results to networks with a fixed 
path-loss from each interferer to the representative receiver, and Section [V] applies the general result to spatially 
distributed networks. Section [VI] contains concluding remarks and a summary of the results presented here. 

The conjugate transpose of a vector a or matrix A are denoted by at and respectively. The determinant is 
represented by |.| and diagonal matrices are represented by A = diag(an, CL22, • • • clnn) where an is the ii-th entry 
of the diagonal matrix A. We use CM(0, v) to represent the circularly symmetric, complex Gaussian distribution 
with mean zero and variance v. The indicator function is denoted by 1_4, which equals one when the condition A 
is true and otherwise. 

II. System Model 

A. General One-to-One Network 

Consider a one-to-one wireless network where there are n + 1 receivers and n + 1 transmitters where each 
transmitter is linked to a single receiver. Let TZ{ denote the i-th receiver and 71 denote the i-th transmitter. Consider 
a representative link between 1Z\ and 71 where % for i = 2, 3, • • • n+1 are co-channel interferers to the representative 
link. 

The representative receiver has N antennas and each transmitting node has K < N antennas. We assume 
frequency-flat fading where the channel between the j-th antenna of transmitting node i and k-th antenna of 
the representative receiver is modeled as ^/Yi9ijk, where 7$ is the path-loss between transmitting node-i and the 
representative receiver and are IID £A/"(0, 1) random variables. We make the standard assumption that nodes 
transmit using Gaussian code -books. 

The representative receiver knows the channel co-efficients between itself and the representative transmitter, and 
also knows the spatial covariance matrix of the interference R = Ki$iK{, where Ki and $1 are defined in 
Section MI- A I Note that receivers can estimate the spatial interference covariance matrix by constructing a sample 
interference covariance matrix from aggregate transmissions of the interferers. 

Each transmitting node knows the channel co-efficients between itself and its target receiver, but not to any other 
nodes. We refer to this assumption as Tx Link CSI. Note that Tx Link CSI can be obtained with low overhead 
in half-duplex systems with reciprocity if channels do not vary rapidly in time since transmitters can use channel 
estimates performed when they acted as receivers in the past, provided that the transmit and receive hardware can 
be accurately characterized. These channel estimates can also be performed by receivers and then fed back to the 
transmitters. 

Additionally, we assume a thermal noise power of No 2 at each antenna of the representative receiver where a 2 
is a constant. The factor N is used to ensure that the limiting SINR is finite so that we can obtain meaningful 
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asymptotic results as N — > oo since the SINR (with a constant noise power) grows at least linearly with N. 
Equivalently, we can assume that the transmit power for each node decays inversely with N as is done in works 
such as [21]. The asymptotic results of this paper apply in the regime as n, N — > oo with the ratio of the number of 
effective interferers (i.e. the total number of independent interfering streams) to receiver antennas nM/N = c > 0. 
All asymptotic results shall refer to this regime. 

B. Spatially Distributed Network 

Consider a circular wireless network of radius R with n wireless transmitters at random IID points in the circle 
such that: 

n = pirR 2 . (1) 

The representative receiver TZi is assumed to be at the center of the circle (defined as the origin) and T\ is an 
additional transmitter at a distance r% from TZi as shown in Figure Q] The n interferers are in links with other 
receivers whose locations do not impact the representative link. Let rj denote the distance between transmitting 
node i and the origin. The path-loss 7, = Gtr^ a with a > 2, which is a standard model for spatially distributed 
networks. 

For the spatially distributed network model, we shall assume a receiver noise power Na 2 where: 

a 2 = a 2 (N~>) (2) 

where a 2 is a constant 

This assumption enables the asymptotic analysis of the SINR in the interference-limited regime as N — > 00. 
Without this assumption, as N — > 00, the optimal receiver will suppress interference to levels comparable to the 
thermal noise which results in the system no longer being interference-limited. Since we focus on the interference- 
limited regime for the spatially distributed network model, a 2 is assumed to be small and does not effect the final 
results significantly. 

III. Parallelized System 
A. Parallelized Transmissions with Link CSI 

Since transmitters do not know the channel coefficients between themselves and receivers other than their target 
receivers, they are not able to encode their transmissions to minimize interference they cause to unintended receivers 
by choice of transmit directions. Since the channels between all pairs of antennas are assumed to be Gaussian 
distributed, no choice of transmit directions is better than any other in terms of interference caused on unintended 
receivers, although the powers allocated to transmit streams can still influence the spectral efficiency of the network. 
Hence, it is optimal for nodes to parallelize the channels between themselves and their targets using an SVD and 
transmit independent data streams on each parallelized channel with some power allocation. We shall assume that all 
transmit nodes parallelize the channels between themselves and their respective receivers and transmit independent 
data streams on M parallel channels where P^ for i = 1, 2, • • ■ n + 1 and j = 1, 2, • • • , M denotes the power 
allocated to the j-th stream by the i-th transmitter. We shall refer to this as the parallelized system. Let the Pij 
be IID over i, i.e. the power allocations of a given transmitter are independent of other transmitters. For a given 
transmitter, there can be arbitrary correlation between the transmit powers it allocates to its M streams. Let the 
Probability-Density -Function (PDF) and CDF of Pij for all i and each j be denoted by fj(x) and Fj(x) respectively, 
and the total transmit power for each node be bounded by Pm, i.e. Y^jLi ^ — ^ M ^ or eac ^ 

Let the N x K matrix «/7jjHy denote the channel matrix between nodes i and j where 7^ is the path-loss 
between 71 and TZj, and Hy is a matrix of IID CM(0, 1) entries (recall that N is the number of receiver antennas 
and K is the number of antennas per transmitter). The spectral efficiency of the representative link is given by (e.g. 
see [26]): 

'Note that the asymptotic SINR we derive for the spatially distributed network is normalized by N a ^ 2 which is accomplished by scaling 
the interference and noise powers by N a ^ 2 . This scaling of the noise requires the noise power to be given by ([2} so that the scaled thermal 
noise value equals Na 2 which is the form of the thermal noise power for the network model of Section Hl-AI 
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Ci = log 2 



n+l 



I + TiHnTiHji Na 2 I + £ 



(3) 



where Tj is the transmit covariance matrix of node-j, i.e., it is the covariance matrix of the signals sent on the 
transmit antennas of node-j. 

Performing an SVD on Hjj yields 



H; 



(4) 



where Ujj and Yij are unitary matrices and Sjj is a matrix containing the squared-singular values of Hjj. Let Xij 
denote the square of the j-th largest singular value of Hji. The spectral efficiency of the representative link can 
then be bounded according to the following Lemma proved in Appendix |A] 

Lemma 1: The spectral efficiency of the representative link can be bounded from above as follows: 



.i 

and below as follows: 



M f -i \ 

Cl < X) lQ g2 ( 1 + JlPlj^Aj ( N ° 21 + Ki*iKj) Uy ) (5) 

.7=1 ^ ' 



Ci>J2 lo S2 ( 1 + TlAjAy-ut (Na 2 I + KjfciKt) u d J (6) 

j=i v / 

where 

*1 = diag(j 2 P21, • • • , 72-FW, 73-P31, • • • ,73-P3M, • • • , 7n+l^(n+l)l> ' ' ' > 7n+l-P(n+l)Af) (7) 

The entries of the N x 11M matrix Ki are IID CN(0, 1) random variables and are defined by (l52l ) in Appendix 
lAl uy are N x 1, unit-norm isotropic random vectors that are mutually orthogonal, and recall that 7^ is the path-loss 
between the i-th transmitter and the representative receiver. For the lower bound, Kj are (N — M + 1) x nM 
matrices with IID CM(0, 1) entries, and iij are (A r -M+l) x 1, unit-norm, isotropic random vectors. 

Note that the upper bound corresponds to the case where the M streams from the representative transmitter do 
not interfere with each other. The lower bound is achieved when the receiver uses M — 1 of its N degrees of 
freedom to completely null the interference from the M — 1 other streams of the representative transmitter. 



B. Asymptotic Spectral Efficiency of the Parallelized System 

The asymptotic spectral efficiency of the parallelized system described in the previous sub-section is characterized 
by the following theorem. Note that Theorem Q] is presented for the general assumptions in Section III-AI and does 
not rely on the spatial distribution of nodes or the path loss model described in Section III-BI 

Theorem 1: Let \£(t) denote the CDF of the path-losses from the interferers to the representative receiver 
72) 73j - "" )7n+i- I n tne limit as n, N — > 00 with the ratio of interferers to receiver antennas nM/N = c, the 
spectral efficiency of link-1 converges with probability 1 to 

M 

C* 1 =Y J ^g 2 (l + X*Pi j 7^) (8) 
i=i 

where (3 is a unique, non-negative solution to the equation: 

Jo l + x/3 
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and H[x) is the limit of the empirical distribution function of the received interference powers given by: 

i M r 

^) = ttE / fj(r)^(x/r)dT (10) 

and Aj is the Umiting value of the j-th largest eigenvalue of the Wishart matrix ^GG^ where G is an N x K 
matrix with IID CAf(0, 1) entries. In particular if N, K — > oo such that K/N = d with < d < 1, then: 

A* = A^ = --- = A M = (1 + ^) 2 (11) 

and if K is a finite constant, 

K = As = • • • = A M = 1. (12) 

Proof: Consider the upper bound for the spectral efficiency of the representative link given in Lemma [U and 
Lemma [2] with m = 0. The upper bound from Lemma Q] can now be written as: 

M 

Ci < J]log 2 (l + SiNRj) (13) 

i=i 

Since 1) SINRj converges with probability 1 to A*Pij7i/3 from Lemma |2l 2) the log function is continuous and 
3) the sum of terms that converge with probability 1 converges with probability 1 to the sum of the limits (e.g. see 
[27]) the RHS of (fT"3l converges with probability 1 to 

M 

J>g 2 (l + A*P li7 i/3)- (14) 
5=± 

Now, let m = M — 1. In this case, the lower bound can be written as: 

M M 

Ci>j2 i°g 2 (i + smR j ) -> E 1 °g2( i + XjPuiiP) ( 15 ) 

j'=i i=i 

The convergence is due to the fact that SINRj converges with probability 1 to A*Pij7i/3 from Lemma[2l independent 
of m. 

Since the upper and lower bounds converge to the same value, we conclude that C\ converges with probability 

1 to Eiiilog 2 (l + A*Pi i7 i/3). ^ ^ ■ 

fi can be interpreted as the limiting Rx array gain SINR as it is the limit of the SINR when the effects of transmit 

power, transmit beamforming, and path-loss are not taken into account. 

Lemma 2: Define SINRj, which can represent either the lower or upper bound of the SINR of the j-th stream 

from the representative transmitter, as follows: 

SINR-, = TiPyAijwt (Vct 2 I + K*iK^ 1 Wj. (16) 

where K is an (N — m)x nM random matrix with IID CN(0, 1) entries with m a non-negative, finite integer, and 
w is an (N — m) x 1 isotropic random vector with unit norm. 

Under the assumptions of Theorem [T] SINRj converges with probability 1 to an asymptotic limit given by: 

SINRj -> jiPijXp (17) 

where the limiting Rx array gain SINR j3, is a unique, non-negative solution for /3 in (© and A* are defined in 

Theorem Q] Note that the limit of SINRj does not depend on m. 

Proof: The proof uses the main results of [23] and is presented in Appendix IB] 

Additionally, we note that the convergence of Ai, A 2 , • • ■ , Am to either 1 or (1 + Vd) 2 is slow and is not a good 
approximation for moderate values of N and K. Recall that d = K/N is the ratio of the number of antenna at the 
transmitters to those at the receivers. Instead, we approximate the j-th largest eigenvalue by the limiting distribution 
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of the eigenvalues evaluated at a fraction N — j of the range between its minimum and maximum values. This 
yields the following approximation: 

Xlj^F^ttN-j + ^/N). (18) 

where Fj 1 {x) is the functional inverse of Fd(x) which is the limiting empirical distribution function (e.d.f.) of the 
eigenvalues of a Wishart matrix ^GGt G here is an N x K matrix of IID CAf(0, 1) entries. Fd(x) is given for 
a general d in Appendix |F] For d = 1, i.e. each transmitter and receiver have the same number of antennas, 

{n+^4x-x 2 +2 arcsin(f-l) 
5- ^ — '- if < x < 4 , 1fY . 
2tT - (19) 
1 otherwise 

Equation ( TTBl can be found by evaluating F7 (a;) numerically. 

IV. Application to Constant Path-loss Systems 
A. Asymptotic Spectral Efficiency 

To test the form of the spectral efficiency described in Theorem [Q we consider two different models for the 
transmit power. In the first model the transmit power used on all M streams by the interferers are constants denoted 
by P, which we call the equal power model. The second model has two classes of nodes and is called the two-class 
model. The first class transmits with power Pi on all streams and the second transmits with power P2 from a single 
stream. The interferers are assigned to classes independently and randomly where the probability that a given node 
is assigned to class one equals q. The two-class model is useful in the context of mixed systems where some 
fraction of transmitters have multiple antennas and the remainder have single antennas. Additionally, the power 
allocated by each user to their transmit streams are correlated, which illustrates the applicability of Theorem Q] in 
systems with correlated transmit powers. 

We assume that the path-losses from all interferers to the representative node equal a constant 7 so that the CDF 
of path-losses ty(x) can expressed in terms of a "step" function with a step at 7, and (© which is the implicit 
equation for the limiting Rx array gain SINR /3 becomes: 

J M ^— ' 1 + xp 



3=1 



For the equal power model, fi(x) = 5(x — P) and (1201) becomes 



*> P + 1 = JgL (21) 

1 + P7/5 



Applying the quadratic formula and selecting the positive term: 



Pep ~ 2a 2 2P 7 + V 4cr 4 + 2P 7 a2 + 4PV ^ 

where we use the notation /3 ep to denote the limiting Rx array gain SINR for the equal-power model. Thus, from 
Theorem 1, the limiting spectral efficiency on the representative link (link 1) is: 

M 

Clep = X) lo g 2 (l + XjP-nPep) (23) 
7=1 

For the two-class model, the marginal PDFs of the transmit powers are: 



and © becomes 



f (x)= \qS(x-Pi) + (l-q)5(x-P 2 ) forj = l 

h[X> \q5(x-Pi) forj = 2,-..,M 



-2a, 1 /3cPi7 , n \ l3cP 2 1 

- (T/3 + 1 = g lTW + (1 " g) M(l + P 2 7/3 ) (25) 
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The exact solution for the limiting Rx array gain SINR for the two-class model is found by solving (125T ). The 
solution denoted by fac, is given by d89l ) in Appendix iGl Hence, the asymptotic spectral efficiency for the two-class 
model is: 

M 

C 12c = lo §2(l + AjPy7lftc). (26) 

B. Monte-Carlo Simulations 

We verified the expressions for the asymptotic spectral efficiency under the constant transmit power and two- 
class models, given by d23l ) and d26l ) respectively, with Ay approximated by (TT8T ). In both cases, we assumed 
a 2 = 1 x 1CT 13 with the thermal noise power equaling Na 2 . 

We simulated systems with the ratio of interferers to receiver antennas n/N = 1 and n/N = 4 with a common 
path loss to the representative receiver of 7 = —125 dB. The representative transmitter had a path loss of —100 
dB to the representative receiver. Each experiment was repeated 1000 times. 

For the constant transmit power model, each transmitting node transmitted P = i W on each of the M data 
streams. For the two-class model, each interferer was class one with probability q = 0.5 where class-one interferers 
transmitted with power Pi = 0.5 W on each of M streams and class-two interferers transmitted P2 = 1 W on a 
single stream. The representative transmitter was always designated as a class-one transmitter. For both models we 
simulated M = 1,2,4 and 8 streams per transmitting node with equal numbers of antennas at all transmitters and 
at the representative receiver, i.e. N = K. 

Figures [2] and [3] show results of the simulations for the constant power model with n/N = 1 and n/N = 4, 
respectively, and the asymptotic spectral efficiency predicted by (|23T ). The points represent a random sampling of 
100 trials from the 1000 trials of the simulation for each N. The standard deviation of the spectral efficiency from 
1000 trials is plotted using dashed lines. The convergence of the spectral efficiency is evident from the figure since 
the points representing different trials of the simulation converge with increasing N. Additionally, note that the 
standard deviation decays with N which indicates convergence in the mean-square sense. For N > 14 antennas 
and 1000 trials, the largest deviation from the asymptotic prediction is less than 15% for n/N = 1 and n/N = 4. 
For N > 25, the largest deviation falls below 10% in both cases. 

Figure [4] show results of the simulations for the two-class model with n = 128 interferers and the asymptotic 
spectral efficiency predicted by d26l ). For N > 14 antennas and from 1000 trials, the largest deviation from the 
asymptotic prediction is less less than 15% for n/N = 4, which is similar to the constant power model, and remains 
below 10% for N > 28. 

V. Application to Spatially Distributed Networks 

A. Asymptotic Spectral Efficiency 

We now apply the results of Section IIII-BI to the spatially distributed network model of Section III-Bi In this 
case it is known from [11] that as the number of receiver antennas N — > 00, the SINR in the interference-limited 
regime for systems without Tx CSI grows as N a l 2 where a is the path-loss exponent. To avoid singularities, we 
define a normalized SINR for the j-th data stream from the representative receiver as follows: 

n jN = N- a/2 SmRj. (27) 

where SINRj is the SINR associated with the j-th data stream from the representative transmitter. 

This normalization is accomplished by scaling the path-loss of the interferers and the thermal noise by N a l 2 . 
The normalized SINR is simply the SINR of this new system with the scaled interference and noise powers. The 
limiting value of the normalized SINR is given by the following theorem which applies Lemma [2] to the scaled 
interference -power model. 

Theorem 2: As the number of interferers n, receiver antennas N, and network radius R — > 00 with nM/N = c 
and n = npR 2 , the normalized SINR for stream-i 77^ approaches an asymptotic limit with probability 1 as follows: 



ViN ^Vi = Pu\*G t r- a (3 



(28) 
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where the limiting Rx array-gain SINR f3 satisfies the following equation: 



2vr 2 p(Gi/3) 




fj(x) X -dx +f3a 2 = l 



(29) 



where b 



a/2 



and E[P^ a ] is the expected value of the transmit power allocated by the interferers to their 

j-th strongest stream, raised to the power -. 
Proof: The proof is presented in Appendix |D] 

In general, d29l ) has to be solved numerically. Moreover, the relationship between parameters such as the number 
of receiver antennas, interferer density, a and the SINR is not clear from d29l . However if we assume that j5 is a 
continuous function of c, a few approximations can be made to yield additional insight into how the various factors 
contribute to the limiting SINR, starting with the following Lemma proved in Appendix HU 



Lemma 3: As b 



for j = 1,2,---M. 



irpN \ 



a/2 



0, if the total transmit power per node is bounded from above by Pm > 0, then 

— — roo 

a —dr / fi(x)x~dx->0. (30) 



From Lemma [3l we note that if n/N is very large i.e. the number of nodes in the network is much larger than 
the number of antennas per receiver, the transmit power limit per node implies that the second term on the LHS 
of d29l ) is small. Furthermore, we assume that the thermal noise power is small, which implies that the third term 
on the LHS of d29l ) is small. Using these approximations we approximate the limiting Rx array gain SINR /3 as 
follows: 



2n 2 p(G t /3)' 



P 



a 
1 



CSC 



2vr\ 
a J 



— ral 



a 



P; 



sm 



2vr\ 
a J 



a/2 



(31) 



Substituting ([3TJ> into ([28) yields a normalized SINR of: 



^ 2 pr\ TIL E 



Pf 



sm 



2?r\ 
a J 



a/2 



(32) 



Since all transmissions use Gaussian code-books, we use the Shannon formula for the spectral efficiency for a 
given link. Writing G a = (^sin(^)) Q/2 and summing the contribution from M streams yields the following 
approximation for the spectral efficiency (with the SINR normalization) of link 1 when N is large: 



M 



I 



§2 



1 + \\iPuG 

a 



V 



1 



2 v^M en 



Pf 



(33) 



where T\ is the length of the representative link and Pu is the transmit power allocated by the representative 
transmitter to its z-th stream. 
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Removing the normalization by N a / 2 we approximate the spectral efficiency of link-1 for large N as: 



M 



( 



Cl ~ l0g 2 



1 + \* U P U G 

a 



\ 



N 



2 t-, 



(34) 



Note that the RHS of (1341 grows with N and hence does not converge. However, the asymptotic spectral efficiency 
of (l34l is a good estimate of the mean spectral efficiency for large N since the deviation of the mean spectral 
efficiency from the asymptotic spectral efficiency is small, as shown in Appendix [EJ Hence, 



( 



£[Ci]«J>g 2 



i=l 



1 + \* U P U G 

c 



\ 



N 



2 v^M pi 



For the equal power model, (1351 ) becomes: 



M 



i=i 



N 



Mirprf 



For the two-class model with Link 1 assigned to class-one, 031 becomes: 



M 



(5 12c «^log 2 l + A^PiG 



i=i 



irprKqMPf + (1 - q)P^ 



(35) 



(36) 



(37) 



and (1351) becomes: 



M 



P[C 12c ] «J]log 2 1 + KiPlGc 



N 



a/2^ 



nprKqMPf + (1 - q)P 2 ~ 



(38) 



B. Monte-Carlo Simulations 

We simulated spatially distributed systems to validate OBI and d38l ). We placed 1000 interferers at random 
locations within a circle of radius selected such that the density of nodes in the network was 10~ 3 nodes per unit 
area and the link-length r\ was such that nprf = 1. The path-loss exponent a was set to 3 or 4 and the thermal 
noise level was constant at 1 x 10~ 13 W. Note that the specific value of the thermal noise power does not play 
a significant role in the interference-limited systems we simulated. In each case the number of antennas at the 
representative receiver N and the transmitting nodes K were equal, and for each N the experiment was repeated 
1000 times. 

For the constant-transmit-power model, each transmitting node transmitted P = i W on each of M data 
streams. For the two-class model, each interferer was class one with probability q = 0.5, where class-one interferers 
transmitted with P\ = 0.5 W on each of M streams, and class-two interferers transmitted P 2 = 1 W on a single 
stream. The representative transmitter was always designated as a class-one transmitter. For both models, we 
simulated M = 1, 2, 4 and 8 streams per transmitting node, with equal numbers of antennas at all transmitters and 
at the representative receiver. 

Figures [5] and [6J show results of simulations of spatially-distributed systems with constant transmit powers, 
a = 4, and SINRs normalized by scaling the interference powers by N a / 2 for 1 and 4 streams per transmitter, 
respectively. The solid lines represent (1331 ) and the dashed lines represent the standard deviation of the simulated 
results. For clarity, only a random sampling of 100 of the 1000 trials are plotted. Figure [7] depicts the simulated 
spectral efficiencies with normalized SINRs for the two-class model and 4 streams per transmitter, with the solid 
line representing (1371 ) with a = 4. In all cases presented here, the spectral efficiencies clearly converge on the 
asymptotic limit as indicated by the distribution of points and the reduction in standard deviation H. 

2 Note that the diminishing standard deviation does not in general imply convergence with probability 1 but it does imply convergence in 
probability. 
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The spectral efficiencies with un-normalized SINR given by (|34l) and (|38l ) do not converge as the spectral 
efficiency increases with N. However, the asymptotic spectral efficiency is a good approximation for the mean 
spectral efficiency for large N as shown in Appendix [E] 

Figures [8] and [9] respectively show the simulated mean spectral efficiencies for the constant power model and 
2-class model with 1, 2, 4, and 8 streams per transmitting node and a = 4. Equations d34l ) and d38l ) are also plotted 
in the figures. 

With constant transmit powers per stream of P/M, Figure [8] indicates that for all simulated cases, the asymptotic 
expression is within 10 % of the simulated mean spectral efficiencies when there are greater than 10 antennas. 
For the two-class model the simulated mean spectral efficiency is within 10 % of the asymptotic expression when 
N > 9 as shown in Figure [9] 

Figure [10] illustrates the simulated and asymptotic mean spectral efficiency with a = 3 and the constant transmit 
power model. For reference, the asymptotic spectral efficiency for a = 4 is plotted using dashed lines. The mean 
spectral efficiency converges to the asymptotic value in a fashion similar to that for a = 4, although the mean 
spectral efficiency is consistently lower for a = 3 compared to a = 4, consistent with (l34l . 

C. Comparison to Systems without Transmit CSI 

When transmitting nodes do not have CSI and transmit with equal power on each stream, the asymptotic mean 
spectral efficiency is known to be [11]: 



where A = Tipr\. This quantity can be interpreted as the average number of interferers closer to a given receiver 
than its target transmitter and was defined as the link rank in [11]. 

To compare the spectral efficiency with and without Tx Link CSI, we plotted the ratio of the asymptotic mean 
spectral efficiency with and without Tx CSI. Figure [TT] shows the ratio of the mean spectral efficiency with Tx Link 
CSI to the mean spectral efficiency without Tx CSI (expressed as a percentage) vs. link-rank A for N = K = 8 
(dashed lines) and N = K = 12 (solid lines) antennas at transmitters and receivers. Note that the gain with Tx-Link 
CSI is highly dependent on A. For instance, for rank-6 links and two transmit streams, Tx-Link CSI doubles the 
spectral efficiency for N = 12. Also note that the increase in mean spectral efficiency can be greater than three-fold 
for high-rank links. 

The increase in spectral efficiency with Tx CSI is greater for large A because the SINR tends to be lower for 
for large A and so the SINR gain provided by the Tx CSI makes a bigger difference inside the log function in the 
expression for the spectral efficiency. The SINR increase in Figure [TT] flattens for large A because at low SINR the 
spectral efficiency of a given stream approaches the SINR because log(l + x) « x for x -C 1. In this regime, the 
ratio of spectral efficiencies with and without Tx CSI does not change significantly with increasing A. The increase 
in spectral efficiency with Tx CSI is lower for a large number of streams due to the fact that the weaker streams 
have SINRs that are close to the SINRs obtained without Tx CSI. 

Figure [TT] thus indicates that a significant (but not an orders of magnitude) increase in spectral efficiency is 
possible with Tx Link CSI that can be acquired with low overhead in duplex systems provided that the transmit 
and receive hardware paths can be accurately characterized. 



A technique is presented for computing the asymptotic spectral efficiency of multi-antenna links in ad-hoc wireless 
networks where transmitters have CSI corresponding to their desired receivers. The transmitters are restricted to 
using M channel modes (which limits the rank of the transmit covariance matrices to M), and the results are 
asymptotic in the regime where the numbers of receiver antennas and interferers go to infinity with a constant ratio. 
The asymptotic predictions are supported by numerical simulations. The simulations indicate that the asymptotic 
expressions are good estimates for the spectral efficiency even when the number of receiver antennas is moderately 
large which is useful since it is possible to place approximately 20 or more antenna elements on a standard laptop 
computer with wave-length separation at a nominal carrier frequency of 6GHz. 




(39) 



VI. Summary and Conclusions 
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In spatially distributed networks the asymptotic spectral efficiency, which approximates the mean spectral effi- 
ciency is found to be dependent on the ratio of interferer density to the number of receiver antennas, as given in 
d34"l ). Thus, as is the case for systems without Tx CSI [11], it is possible to maintain a constant mean spectral 
efficiency per link if the number of antennas per receiver is increased linearly with interferer density. 

Additionally, we found that the spectral efficiency in the network can be increased in certain cases if each 
transmitter transmits fewer data streams unlike MIMO links in AWGN channels where to maximize capacity, 
nodes should transmit data on all their channel modes with a water-filling power allocation. Figure [8] illustrates 
this potential advantage of transmitting fewer streams when the number of antennas is relatively small; a similar 
observation was made in [11] for systems without Tx CSI. 

Compared to systems without Tx CSI, we find that the asymptotic spectral efficiency with Tx-Link CSI can be 
several times larger, where the benefit of using Tx-Link CSI increases for longer links, denser networks, or both, as 
illustrated by Figure [TT] For instance, with four transmit streams and 12 antennas at transmitters and receivers, Tx 
CSI can double the spectral efficiency when the link-rank A = 3, where A = irprf, which captures link length and 
interferer distribution. Since Tx-Link CSI can be estimated in duplex systems with reciprocity without a significant 
increase in overhead, Tx Link CSI can be useful, particularly in dense networks with long links. 
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Appendix 



A. Bounds on the Spectral Efficiency of Parallelized System 

The spectral efficiency of link i between 1Z{ and % is given by (see e.g. [26]): 



Ci = log 2 



n+1 



I + 7 «H«TiHi ( iVcT 2 I + Yl 'U u .i> T .i H ) 



(40) 



where 7^ is the path loss between transmitter T% and TZj, and Tj is the transmit covariance matrix of Tj, i.e., it 
is the covariance matrix of the signals sent on the transmit antennas of node j. Recall that ^/j^jHij is the N x K 
matrix of channel coefficients between the antennas of % and TZj. 

Since transmitters do not know the channels between themselves and unintended targets, they cannot choose their 
transmit covariance matrices so as to transmit in spatial directions that reduce interference to undesired receivers. 
Hence, they should transmit in spatial directions that maximize the data rate on their individual links. Note that 
their choice of transmit powers to allocate to their streams can still influence the spectral efficiency of other links. 

It is known that without knowledge of the quantity in the parenthesis in (l40l . to maximize the RHS of (1401 the 
i-th transmitter should use the following transmit covariance matrix: 



with 



V,P,V 



diag(Pa,JW- - ,PiM,0, 



(41) 



(42) 



where Pj^ is the power allocated to the A;-th stream by the j'-th transmitter. V, is a unitary matrix defined by taking 
the SVD of Hjj such that = UjXljVt. Sj here is a diagonal matrix of the singular values of Hj with A^- equal 
to the square of the j-th largest singular value, and Uj is another unitary matrix. 

Assuming that all transmitters use the same strategy, substituting (|4~TI ) into (1401 yields: 



d = log 2 



n+1 



I + 7ii H«ViP i v{Hi [ Na 2 I + 7<jHi<VjPjVtHt 



(43) 
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Note that random matrices with Gaussian distributed entries maintain their statistical properties when multiplied 
by unitary matrices. Thus, we can write (|43T ) as 

-l 



Ci = lop 



n+l 



I + 7 «H«V i P f V{l4 ( Na 2 I + £ Tij'HjiPjH 



where Hjj are distributed identically to Hy. 

Substituting Hii = UiEjVj and noting that V;Vj = I: 



Ci = loo 



log 2 



n+l 



I + 7ii U i S i P i E}ut [ Na 2 I + TiiHjiPj-H 



n+l 



t 

-1 



I + 7«S i P i S,tut [ iVa 2 I+ £ TijHjjPjHt. J u, 



(44) 



(45) 



The steps from (|4Q| > to (1451 ) are standard and can be found in [26]. 

1 ) Upper Bound on Spectral Efficiency: Consider the spectral efficiency of the representative link, link 1 . For 
notational convenience, define: 



Q = Uj \^Na 2 I + £ 7i J H i iP,HU Ui 
With qjk denoting the jk-th entry of Q, we have: 



(46) 



(47) 



where ui^ is the A;-th column of Ui. 
We can write (I43T) as: 



Ci = log 



I + 7ll £ 1 Pi£ t 1 Q 



(48) 



Si contains the singular values of the channel matrix of the representative link with Ay equal to the square of 



the j-th largest singular value. Thus, from (1421) the j-th diagonal entry of I + 7nEiPi53|Q is 1 + JixXijPijqjj 



(49) 



Hence, by the Hadamard inequality (see e.g. [28]) we can bound (1481 as follows: 

(K \ M 

Yl (1 + ixiXxjPijqjj) = Y log2 ( x + " i '■ 'V/ l'<j'i.jj • 

where the last step uses the fact that Py = for j > M. Substituting (1471 yields: 



d < £ log 2 1 + -:iAo/', ; u |; Na 2 I + £ " :,H ; , I', ] u y 

J'=l \ \ J=2 



Using some matrix manipulations, (1501 ) can be written as: 



M 



i=i 



Cx < J2 lo S2 1 + - 1 ; A , , , ii ( / (iVa 2 I + Kx&xKU u y 



(50) 



(51) 



where <&i is given in ([7]) and 



H 21 H 31 



(n+l)l 



(52) 
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where the N x M matrix comprises the first M columns of Hji. 

Note that 3>i is a diagonal matrix containing the received powers from each stream transmitted by each interferer 
and Ki is a matrix of IID CN(0, 1) entries. 

2) Lower Bound on Spectral Efficiency: Since the transmit covariance matrix of the j-th transmitter is VjPjVj, 
the vector of received samples at the antennas of the representative receiver at a given sampling time can be written 
as follows: 



n+1 

y = HnViPfxi + J] N -— H,-. Pjx; + n 

3=2 



(53) 



where Xj contains unit variance samples from Tj in its first M entries with the remainder being zeros, and n is an 
N x 1 vector of IID CjV(0, a 2 N) entries representing the thermal noise. 

We shall find the lower bound by using a specific, suboptimal procedure where the representative receiver decodes 
each stream from the representative transmitter individually. When decoding the i-th stream from the representative 
transmitter, the receiver will use a fraction of its degrees of freedom to completely null the interference from the 
other data streams of the representative transmitter. 



Suppose that the receiver multiples y by Uj, which yields: 



n+1 



y = Ujy = IjjHuViPf Xl + uj £ v /7iJH il P|x j + ujn 



3=2 



n+1 



UjUiSiV|ViP|xi + £ v^i-ulH^Pjx^ + Ujn 



3=2 



n+1 



EiPf X! + £ v~H ; ;l';x ; + n 



3=2 



n+1 



*i + v/TliHjiPjxj +n 

3=2 



(54) 



where H^i = Tj|Hji is a matrix with IID CJ\f(0, 1) and n = Ujn contains IID CjV(0,a 2 N) entries. Note that 



/ VAnPiixn \ 



xi = EiPfxi 







(55) 



V o / 

where the zero entries appear because P\j = for j > M and x\j is the jth data sample of the representative 
transmitter, i.e., the jth entry of the vector xi. Note from d54l and (|55T ) that the data samples from the representative 
transmitter are parallelized in y so each sample of y contains information from a single stream of the representative 
transmitter. 

To decode the ith data sample from the representative transmitter the receiver constructs a vector yj by discarding 
samples that contain information from the remaining M — 1 streams of the representative transmitter. The (N — 
M + 1) x 1 vector yj is defined as follows: 

/ m \ 

V(M+l) 



y; 



(56) 



V VN J 
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where yj is the jth entry of the vector y. Thus, yj can be written as follows: 



\ 







n+l 



+ X! \ / ~ H " ,, j x / + n 

i=2 



(57) 



where Hjj is an (N — M + 1) x K matrix of IID CN(0, 1), which equals the matrix Hji with the rows 
1, • • • , i — 1, i + 1, • • ■ , M removed. The (iV - M + 1) x 1 vector n contains IID CW(0, <r 2 iV) noise samples. 

Suppose now that the receiver picks an (N — M + 1) x (N — M + 1) matrix U with uniform probability from 
the group of all unitary matrices and multiples it with yj. This yields the following: 



n+l 



fi = Uyj = y/ XuPuiiiXi + ^2 Xj + n. 

J =2 



(58) 



where Uj is the ith column of U. Since U is unitary, EL-^ = UHj! is still an (N — M + 1) x K matrix of IID 
CAf(0, 1) entries and n = Un is an (N - M + 1) x 1 vector of IID CW(0, a 2 N) noise samples. 

The receiver then uses a linear MMSE receiver on the vector yj to detect for which the SINR is well known 
(e.g. see [7]): 

(n+l \ _1 _ x 



where 



K, 



H 2i H 3i ' ' ' H (n+l)i 



(60) 



where is the matrix comprising the first M columns of Hjj and $i is defined in ©. 

Since all transmissions use Gaussian codebooks, each stream can support a spectral efficiency of 



log 2 + XuPuu] (K^iKt + iV,j 2 I 
Summing the contributions of M streams gives the lower bound. 



u 



(61) 



B. Proof of Lemma |2] on the Limiting SINR per stream. 

Note that in (fT6l) is an isotropic random vector with unit norm and (e.g. see [22], [29], [30] or [31] Appendix 
A.l.) can be expressed as[l : 

Wj = rp-TTgj (62) 
where the entries of gj are distributed as IID CM(0, 1), which means that (fT6l ) can be written as: 

smRj = ^# 7iPi ^ Ai ^ Sjt (^ 2l + ^ Ki * iK 1 gj m 

If the e.d.f of the entries of <l?i converges with probability 1 to a limiting function H(t), the factor 

+ ^ Kl * lK (64) 

was shown by [23] to converge with probability 1 to an asymptotic limit as n — > oo. The following lemma proved 
in Appendix O shows that the e.d.f. of the entries of 3>i indeed converges with probability 1 to an asymptotic limit. 

3 Note that expressing isotropic vectors in this manner was used in [22] in their analysis of the fluctuation about the mean of the SINR of 
random CDMA systems with signature vectors comprising IID entries. 
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Lemma 4: Let H n {x) denote the e.d.f of the interference powers for any given n, i.e. H n (x) is the e.d.f of the 
entries of the diagonal matrix Then as n — > oo in the manner of Theorem [TJ 

i - r 

H n {x) -> H{x) = JjJ2 fi( T M x M dr (65) 
i=i 17 

JL n-t ( =.2t i TT_ rT» irt 



From Lemma @] and the main result of [23], the term -^gt ( a 2 I + Ki$iK|J gj converges with probability 1 

to an asymptotic limit which we define as /3. From [21] igt (a 2 I + Ki$iK|^ g converges in probability to (5 
which is a unique solution for (3(z) in the equation: 

when z = —a 2 . Since [23] proves convergence with probability 1 to an asymptotic limit as n — > oo, and [21] 
proves convergence in probability to the solution to d66l ), we conclude that ^gt [a 2 I + Ki<I>iK{^ gj converges 
with probability 1 to the solution to 
Additionally, note that 



jvll&ll 



with probability 1 by the strong law of large numbers. 

If K — > oo as N — > oo with K/N = d, then 4° Ay for j = 1, 2, ■ ■ ■ ,M is known to converge with probability 1 
to an asymptotic limit A* = (1 + Vd) 2 for j = 1, 2, • • • , M (e.g. see [32]). If K is a finite constant, by the strong 
law of large numbers -^Xij for j = 1, 2, • • • , M converges with probability 1 to unity. 

Hence, each random term in (fT6l ) converges with probability 1 to a non-random value. It is known that for 
sequences of random variables X n and Y n , if X n — > X and Y n — > Y with probability 1, then X n Y n — y XY with 
probability 1 as well (e.g. see Theorem 5.21 in [27]). Hence, (fT6l) converges with probability 1 to 



7 iPi;A*/3 (67) 

C. Proof of Lemma on the convergence of the Interference Powers 

Let pij = jiPij, i.e. pij is the product of path-loss to the representative receiver and the transmit power on the 
j-th stream of the i-th transmitting node. Recall that N, n — > oo with n/N = c/M > 0. Thus we have the limiting 
e.d.f. of the received interference powers: 

. M n+l 

H(x) = lim H n (x) = lim — V - V lr p <a;} (68) 



M ^ n 

j=l i=2 



1 M 

m22, E ^{ph^}] with probability 1 



M 



1 M 

^ Pr{py < x] 



M 



(69) 
(70) 



where the step from d68l ) to (|69l follows from the strong-law-of-large-numbers. 
Since p^j are distributed identically for all i, 

Pr{ PlJ < x} = J fj(r)Pr{ PlJ < z|i^ = r}dr (71) 

Substituting p^j = jiPij yields: 

Pr{pij < x} = J / j (T)Pr{7 i r < = r}dr 



(72) 



fj{T)9(x/T)dT (73) 
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Thus we have: 

H(x) = lim H n (x) = [ /j-(r)*(x/r)dT (74) 

D. Proof of Theorem [2] on the limiting SINR for the spatially distributed network model 

Recall that the path losses from each interferer and the thermal noise are scaled by TV"/ 2 for this model. The 
SINR on the j-th stream from the representative transmitter can be bounded as follows: 

jiPijXiju] (Na 2 I + KjSiK})" 1 uj < vnj < -;/V\| ; u |; (V<r 2 I + Ki^Kj)" 1 u y (75) 

where u.j, Uj, Kj and Ki are as defined in Lemma [T] Recall that the upper bound is attained if the M streams 
from the representative transmitter do not interfere with each other and the lower bound is attained by perfectly 
nulling the interference of streams 1, • • • , (i — 1), (i + 1), • • ■ , M from the representative transmitter. 

If this new system with the scaled path-losses and noise power meets the conditions of Lemma |2l then the upper 
and lower bounds in d751 ) will converge to the same asymptotic limit, implying that the normalized SINR rjNj 
converges to that limit as well. The rest of this section is devoted to showing that the system model meets the 
requirements of Lemma |2] and to finding the limiting value of the upper and lower bounds in d75l ). 

We start by showing that the empirical distribution of received interference powers converges with probability 1 
to an asymptotic limit H{r). For this network model and a given n, N, and R, let the CDF of the path-losses be 
denoted by ^at(x). Note that the interference power is scaled by N a / 2 here. 

V N (x) = Fr{N a / 2 G t r7 a < x } 
Pr < r„; > 

R 



G t 

xN-<*/ 2 * 
G t 



Substituting ©, c = nM/N, and b = 



R 2 {o<(^iyKn} 

a/2 



(76) 



n ( xN~ a/2 



7rp V G t 



77 p 



Nirp i 


' X 


n ' 


\G t/ 


irpM 


(- 








L {G t fe<x<oo} 



^(x) 



Hence ^^(x) is independent of N (although it depends on n/N which is a constant) and Lemma [2] holds for the 
spatially distributed model with scaled interference powers. The remaining steps are to find H(x) and its derivative, 
which are used to evaluate (fl6l ). which gives the limiting normalized SINR. 
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3=1 

Jj-f E / M T ) ( X " ^ (^f) " 1 {G t 6<x/r<oo} > j dT 



3=1 







Ef bG t 2 

— M — M 

* pG? e r ^ + ^ e r ^ * 

-•_ 1 JO CI" • , Jt^- 



The derivative of is: 



ex. J=1 Jo ex. i=1 -_ 

^E^ + ^E/ /<w^ 



drc acx + . i=1 acx + « i=1 7^- 
Substituting (1771 ) into (O and scaling the thermal noise power by N a / 2 such that a 2 = a 2 yields: 



A M _ „ ^4 AT „oo 



2o , 1 a f x SfpGff p . 2, 2vrpGf ^ f°° , , , 
<r 2 (3 + l = Pc — — 1+ 2 > E[P a ] ^— V > / /; t t 



^^rAL^M^r^Lr^^ (78) 

« c Jo 1 + « c ~-f Jo 1 + y^- 

J — 1 J— 1 bG t 



From the proof of Lemma 1 in [llfl 



dx = P«- l itcsc ( — J (79) 







1 + x/3 \ a / 



Substituting 479]) into (1781 ): 



ac v . 7 .. 



a / ac Jo l + x/3 

Rearranging terms yields d29l ) which completes the proof. 



Note that the lower limit of the integral in Lemma 1 of [11] is greater than zero, however its proof clearly allows the lower limit to be 
zero 
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E. Approximating the Mean Spectral Efficiency by the Asymptotic Spectral Efficiency 

Consider the deviation of the mean spectral efficiency from the asymptotic spectral efficiency on the i-th stream 
If we assume that E [log 2 (rim)] is bounded for all N, then by the bounded-convergence theorem (see e.g. [27]), 

'iV~f + rj Ni 



\E [log 2 (1 + N* rim)] ~ log 2 (1 + N'th) | 



E 



log 2 



NS +T]i 



0. 



(80) 



This implies that the deviation of the mean spectral efficiency from its asymptotic value decays to zero, i.e.: 

M M 

J2E[log 2 {l + Nf Vm )] -^log 2 (l + 7Vt^) ^0. (81) 
i=i i=i 

Thus for large N, the mean spectral efficiency is well approximated by the asymptotic spectral efficiency. 
To show that E [log 2 (rjm)] is bounded for all N, consider: 

POO 

|#pog 2 (i7jvi)]| <-E[|log 2 (r/jvi)|] = / Pr {|log 2 (vm)\ > %} dx 

Jo 

POO POO 

= / Pr {vm > 2 X } dx + / Pr { VNi < 2~ x } dx 
Jo Jo 

is bounded because rjm is bounded from above. Now, consider the second 



The first term on the RHS of 
term on the RHS: 



Pr{ VNi <2- x } <PT{fj Ni <2- x } 



(82) 



where r)Ni is the normalized SINR of the representative link in an infinite wireless network with N antennas at 
the representative receiver, no CSI at the transmitters, and any suboptimal linear receiver. 

Consider a suboptimal receiver which uses only one of its antennas selected at random if N < for some 

< 9 < i. In this case from [33], 



Pr {fjm <2~ x } = 1- exp i-G 1 N2~ x - G 2 N2~ 
< 1 - exp (- (NGi + NG 2 ) 2"' 



(83) 

where G\ and G 2 are positive parameters independent of x. The integral of ([831 w.r.t. x is finite. 

If iV > iff], the receiver uses the partial-zero-forcing algorithm of [16] with k = 9N degrees of freedom used 
for zero-forcing. In this case, we have from [16]: 



Pr {fj Ni <2~ x } < 



N~2 2- x (G 3 (ON - [f I) 1 ' + G 4 iV 1 - Q/2 a 2 ) 



(l-e)N-i 

l 2~ x G 3 (ON- [fl) 1 ~ f +G 4 a 2 



< 



(l-0)-l/N 
'H- X G 3 (ON - [f]) 1 "^ + G 4 a 2 



Ni 



1 + TIT 



2~ x Gi 9 



2~ X G 3 9 



< 



i + T 4- 

I 2 

1 ra 

rirtn 1 2 



+ G 4 a 2 



1 



1 + 



(84) 



where G3 , G4 are positive terms independent of N and x. 

The integral of the RHS of (l84l is clearly finite which implies that E [|log 2 (??jVi)|] is finite confirming (I8TT ). 
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F. Asymptotic Empirical Distribution Function of the Eigenvalues of Wishart Matrices 
It is known that (e.g. see [32]): 

dF fy^ = max (o> 1 - d )${y) - («2 - y){y - ax)i {ai<y<a2} (85) 

where a± = (1 — \fd) 2 and a 2 = (1 + \fd) 2 . Taking the integral with respect to y from to x, we find that for 

ai < y < a 2 



1 1 1 /t r- r 1 /ai + 02 — 2y\ 

Fd{y) = o( a i + °2) - 7Vaia2 + — V ( a 2 - y)(y - ax) + -7-\fl\ + a 2 ) arcsm 

O 4 27T 47T V ai — a 2 / 

1 / 2 a\a 2 - a 2 y - aiy 

— yja\a 2 arctan — 

2vr V 2^0102(02 - y)(y - ax) 



1 , , / 2aia 2 - a 2 y - axy \ 

+ 7— Vaxa 2 arctan — = (86) 



and 



F d (x) = h axi0A ~ d) if0 ^«i (87) 
I 1 if a 2 < x 

G. Asymptotic Spectral Efficiency for Constant Path-loss and the Two-Class Model 
For simplicity, we express (|25T ) as 

H/3 3 + T 2 /3 2 + T 3 /3 - 1 = (88) 

where 

Tx = a 2 PxP 2l 2 

(1 - q)cPxP 2 ~f 2 2 2 2 

T 2 = — h cqPxP 2 T + {? ~ P\l)P2l + cr -P17 

™ d 0.-2, (i-g)ciVy . D 

?3 = -P27 - + o- H 1- QcPxl 

Additionally, we define the following to simplify notation: 

T 4 = 2T 2 3 - 9TxT 2 T 3 - 27T 2 



T 5 = ^T 2 - 4(T 2 2 - 3TiT 3 ) 3 
The solution to (l25l) is found using the standard formula for the roots of a cubic polynomial, which yields: 

/32c --3T\ +_ 6T 1 — {-^2— ) + -&T\—{-^2—) (89) 
//. Proof of Lemma\3\ Used to Simplify Limiting Mean Spectral Efficiency 

Note that /3 must mono tonic ally decrease with b = {^—^j , and that the maximum transmit power per stream 
is Pm since the total transmit power is bounded by Pm- Thus, fi(x) = for x > Pm and so for any i, 

fi(x)x~dx dr 

fi(x)x^dxdT (90) 







1 + r/3 

G t b 

bPu T - 







l + r/3 



Gtb 



bP M T - 



*/ TT70 P " dT (92) 

which goes to zero as b — > 0. The step from (l90b to (|9TT) is because the maximum transmit power is Pm and (l92l) 
follows because PDFs integrate to unity. 
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■ Representative Receiver 
• Representative Transmitter 
- Representative Link 




Fig. 1. Illustration of wireless network with representative link. 
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Fig. 2. Simulated and asymptotic spectral efficiency vs. number of antennas with constant path-loss, constant transmitter powers, and 
n/N — 1. The dashed lines represent the standard deviation of the simulated spectral efficiencies. 
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Fig. 3. Simulated and asymptotic spectral efficiency vs. number of antennas with constant path-loss, constant transmitter powers, and 
n/N — 4. The dashed lines represent the standard deviation of the simulated spectral efficiencies. 
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Fig. 4. Simulated and asymptotic spectral efficiency vs. number of antennas with constant path-loss, transmitter powers from the two-class 
model, and n/N — 4. The dashed lines represent the standard deviation of the simulated spectral efficiencies. 
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Fig. 5. Simulated and asymptotic normalized spectral efficiency vs. number of antennas for spatially distributed networks with constant 
transmit powers and one stream with a = 4. 
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Fig. 6. Simulated and asymptotic normalized spectral efficiency vs. number of antennas for spatially distributed networks with constant 
transmit powers and 4 streams with a = 4. 
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Fig. 7. Simulated and asymptotic normalized spectral efficiency vs. number of antennas for spatially distributed networks with the two-class 
model and one transmit stream with a = 4. 
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Fig. 8. Simulated and asymptotic mean spectral efficiency vs. number of antennas for spatially distributed networks with constant transmit 
powers and a = 4 . 
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Fig. 9. Asymptotic and simulated mean spectral efficiency vs. number of antennas for spatially distributed networks with the 2-class model 
and a = 4 . 




Fig. 10. Simulated and asymptotic mean spectral efficiency vs. number of antennas for spatially distributed networks with constant transmit 
powers and a — 3. The asymptotic spectral efficiency for a = 4 is shown by the dashed lines. 
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Fig. 11. Percentage increase in asymptotic mean spectral efficiency of systems with Tx Link CSI over systems without Tx CSI vs. -Kpr\ 
for N = 12 (solid lines) and N — 8 (dashed lines) antennas per node. 



